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X-ray Navigation - XNA V 


• Millisecond Pulsar X-ray Navigation 

- Pulsars are neutron stars that appear to pulsate across the 
electromagnetic spectrum, but are bright enough in X-rays to 
enable XNAV instrumentation & solutions 

- Some millisecond pulsars (MSPs) rival atomic clocks in long- 
term timing stability 

- Pulse phase or time of arrival deviations can update spacecraft 
dynamics model along line of sight 

- Observing multiple MSPs can yield 3D orbit data 

- GPS-like positioning & timing throughout solar system & 
beyond 
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XNAV Development History 



XNAV has rich history beginning with discovery of first radio 
pulsar 

- Significant body of published research 
Naval Research Laboratory (NRL) (1999-2000) 

- Unconventional Stellar Aspect (USA) Experiment 
DARPA XNAV Project (2005-2006) 

- Ball Aerospace collaborated with Microcosm Inc. 

- Algorithms, Infrastructure 

- Detector and Pulsar modeling studies (NRL) 

- Modulated X-ray Source (MXS) developed, Gendreau 
DARPA XTIM (2009-2012) continuation DARPA XNAV, led by 
Lockheed with Ball Aerospace 

- Used Large Area Collimated Detector 

NASASBIRs with Microcosm LOCKHEED MARTIN 

NICER / SEXTANT selection 04/2013 

- SEXTANT team deeply involved in prior programs 

- Evolution of XNAV detector ideas shows NICER XTI 

(concentrating optics/ silicon det) to be practically ideal *..mm * 






s \^!Ty 


Prior work has set the stage for SEXTANT to perform the full 
onboard XNAV OD 



Microcosm 
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NICER/SEXTANT Companion Mission 



NICER - Neutron Star Interior Composition ExploreR 

- Fundamental investigation of ultra-dense 
matter: structure, dynamics, & energetics 

- Launch in Oct 2016 on Space-X Dragon 

- 18 Month mission on Express Logistics Carrier 
(ELC) 

- X-ray Timing Instrument (XTI) 

• X-ray (0.2-12 keV) concentrator (single-bounce) optics 
and silicon-drift detectors 

• High precision time tagging (300 ns absolute) 


• Large effective area (>1800 cm 2 ) 


• X-ray detectors with high quantum efficiency and 
spectral resolution 

SEXTANT -XNAV Flight Demo 

- STMD/GCDP funded technology enhancement 

- 1 st demo of real-time, on-board pulsar XNAV 

- Leverage NICER Hardware with enhanced flight 
XNAV software 


X-Rays 



28V 


RS422 Link for 
Commands & data 


GPS Pulse 
Per Second 


Block diagram of NICER XTI 





XNA V Demonstration Concept 





NICER 

Instrument 


X-ray photons 


GPS XNAV 
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Folded Light Curve for PSR-J0534+2200 



1 folded counts 
-A- fit pulse 





SEXTANT 

algorithms 
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• SEXTANT Primary Objective 

Provide first demonstration of real-time, on-board X-ray Pulsar Navigation 

— Primary Objective: 10 km orbit determination accuracy, worst direction, 2 weeks 
— Implement a fully functional XNAV system in a challenging ISS/LEO orbit; NICER compatible 

- Advance XNAV technologies 

• Secondary Objectives 

- Repeat XNAV demo with new NICER data & new MSP discoveries 

- Validate & enhance the unique Goddard XNAV Lab Testbed (GXLT) 

- Use SEXTANT data & GXLT to study real-world XNAV scenarios 

- Evaluate candidate photon processing & navigation algorithms & develop new techniques 

- Study utility of pulsars for time keeping & clock synchronization 

• Planned Experiments 

- 2-week period observing 3-5 pulsars early in the mission (primary) 

- Opportunistic experiments 

- Ground experiments using collected photon data 

• Stretch Objective 

- 1 km orbit determination accuracy, worst direction, using up to 4 weeks of observations. 
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SEXTANT System Architecture 


NICER XTI 



• 1 . 

t ti tint 


GPS time stamped 
Photon + Background 
events, GPS telemetry 


GPS time stamped 
Photon + Background 
events 


Ground System 
Software 


Ground System 
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Algorithms 
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Pulsar almanac 
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Pulsar almanac entries (timing models, 
templates) and alert messages 
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Test stimulus 


Pulsar data from observatories 
(Arecibo, FERMI, etc.) 


Testbed 


Goddard XNAV 
Lab Testbed 
(GXLT) 


Proving ground for algorithm 
development, validation, test 


Test stimulus 
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NICER X-ray Timing Instrument ( XT I) 


56 co-aligned X-ray concentrator 
optics and associated Silicon Drift 
Detectors (SDDs) in Focal Plane 
Modules (FPMs) 

7 Measurement/Power Units 
(MPUs) 

The FPMs detect X-rays arriving 
from the concentrators 



antennas (2) 


0.93m 


itrator Plate 
Star Tracker 


NE2227 


Electronics 
Plate 


Electronics 
(MPUs, MEB, GEB, Star tracker DPU) 


Paraffin Panels 


Detector 


Radiator 


MPUs time-tag and packetize 
photon events 

< 300 nsec absolute time resolution 

> 1800cm 2 effective area 

Moderate (CCD-like) energy 
resolution 




Stowed deploying , and observing 
at ELC mounted location on ISS 



• XNAV Flight Software (XFSW) Sequence Flow: 

- XTI detects events from sequential pulsar observations, output via MPU 

- Pre-process filtering of photon events and buffer collection to obtain 


statistically sufficient observations of a single pulsar 

- Batch process events to extract single measurement of phase, Doppler, count 
rates 

— Navigation algorithm (GEONS EKF) blends models of S/C dynamics with phase 
& frequency measurements to update spacecraft state estimate 


• Ground system maintains pulsar almanac used by XFSW 









Photon Event Arrival Model 


Crab lightcuve and assoc, photon arrivals 


The arrival times of X-ray photons at the 
detector are modeled as a non- 
homogeneous Poisson process (NHPP) with 
quasi-periodic rate function: 

A(t) = A (<j>(t)) = /3 + ah(cj>(t)) 

• h is the known pulsar lightcurve 

• periodic in [0,1) integrating to 1 with 
minimum value 0 

• b is total background rate 

• a is the mean signal rate 

The probability of k events arriving in the 
interval (t 0 , t b ) is given by 

P (M„. et)) . - «(..))■ 


k\ 


A(t) = /o A (s)cfe 



Arrival times follow a NHPP 


Folded Light Curve for Crab pulsar 


■folded photon events 
— reconstructed pulse 
—zero phase reference pulse 



0.015 0.02 

Time (sec) 


Lightcurve recovery by "pulse folding" 
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(j){t) = 0 o (t-r(t)) 


For a Geocenter RefObs: 


- 0 o is the phase evolution at a reference 
observatory (RefObs) 

- provided by the standard pulsar timing software 
TEMP02 (see Hobbs et. al v Edwards et. al.) 

• Extremely high fidelity models 

• Provides convenient piecewise polynomial 
approximations to the full timing model 4> 0 

- r(t) the prop time of the pulse wavefront moving 
from the detector to the RefObs at speed c 

• If the RefObs is close to the detector 
(Geocenter) T ^ _ *(0 ' n 

c 

where x(t) is the detector coordinates in a 
frame centered at the RefObs and n is the 
direction to the given pulsar 

- If the RefObs is not so close to the detector (e.g., 
SSB), parallax and Solar Shapiro delay terms are 
needed as well 

- SEXTANTs algorithms are based on a Geocentric 
RefObs but can support SSB RefObs 


0(0 = 00 (t 


x(t) • n 




~ 0o t 


c 

~x(t) ■ n 


+ 4*0 \ t 


(x(t) from filter) 
x(t) ■ n\ Sx(t ) • n 


=: 0(t) + e(t) 

where 

Sx(t) = x(t)-x(t) 

For SEXTANT, we assume 

e(t) ~ q+ f(t - f a ) (£? and/are 

constants) 


Alternate models are possible but this 
linear model works well in simulations 
(with short observation intervals) 
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Pulse Phase Estimation Algorithm 


1. Observe arrival times {7"/ f }^ =1 during a fixed interval [t a , t^\. 

A 

2. Determine estimates (q, f) of the parameters (q, f) in the 
model e(t) = q + f(t — t a ) by maximizing the log-likelihood 
function with {7"/ c }^ =1 

N 

L(§) = J2'og\(4>m + q+ f( T k - t a j) - (a + P)(t b - t a ) 

k=l 

with 6 = (q, f, a, /3). 

3. Form phase and frequency estimates 

A ■ 

4>(t) = 4>(t) + q + f(t — t a ) and (/>(t) = </>(t) + f , respectively. 


13 



Phase Estimate Accuracy 


The Cramer-Rao lower bound (CRLB) gives limit of achievable accuracy of any estimator with fully specified 
statistical model. 

We can compute CRLB for NHPP model assuming a, (3 known 

- For XNAV the CRLB is smallest for short period, for highly peaked pulsars with high SNR 
ML estimators asymptotically (many photons) achieve CRLB 

CRLB useful for 

- performance estimation 

- determining target observation times 

There are other error sources including those associated with model mismatch 

- best you can do with a statistical model 
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GEO NS Navigation Filter 


• Goddard Enhanced Onboard Navigation System (GEONS) 

- Flight-proven, Award-winning, NPR 7150.2 Compliant UD-Factorized Extended 
Kalman Filter 

- Estimate multiple spacecraft absolute and/or relative states 

- Enables data fusion and regime independence 

- Earth, Moon, LPOs, Deep Space 

- GPS, TDRSS, DSN/USN/GN, Crosslink, Celestial Object, Accelerometer, and XNAV 
measurements. 

- Used (and will be used) on Terra, EO-1, GPM, MMS, NICER/SEXTANT 

- Licensed to Orbital, Ball, ITT, Moog/Broad Reach 

• SEXTANT team added an XNAV measurement model to GEONS 



where v is noise 


connects pulse phase observable with spacecraft state 
• Differentiate to get a frequency measurement equation 
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Level 0 simulation Level 2 simulation 


• Software only XNAV measurement simulation • 

• Useful for long term studies (deep space • 

trajectories, etc) • 

Level 1 simulation 

• Software only photon event simulation • 

• Photon Processing algorithm implemented for 

measurement generation • 

• Primary mode of development for SEXTANT 


Hardware-in-the-loop simulation 
Test-as-you-fJy 

Use the Modulated X-ray source (MXS) to 
generate the photon events 

X-ray detector & electronics time-tag the 
photon events 

Useful for testing flight hardware 



Legend 
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I/O Data 

j Alternate simulation paths 

Software accelerated mode: 



Software 

1 measurement-level simulation 



Process 

' Hardware mode (real-time only) 



Hardware 

Component 

Software accelerated mode: 
photon-level simulation 




> 



Measurement 
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• Achieving performance 
levels that meet 
SEXTANT requirement 
(10km worst direction 
after 14 days) with 
margin. 

• Three pulsars: 
B1937+21,B1821-24, 
J0218+4232 

• Large initial errors added 
to the state 

• Rapid divergence occurs 
without processing XNAV 
measurements 


GEONS filtered pos/vel errors 

^ 15 r 

E 
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Generate and maintain the pulsar 
almanac 

- Timing models 

- Profile templates (light-curves) 

- NICER data after sufficient data 
collected 

Maintain and update application 

- Pulsar upload tables 


Ground Radio & Space 
Observatories 



i 


Raw Timing 
Data Analysis 


Generation of 
Timing Models, 
Templates, Rates, Etc. 


SEXTANT Ground System 


• Monitor performance 

— Trending 
- Alerts 

• Driven initially with radio 
observations 


SEXTANT Pulsar 
Catalog 




Products for 
XFSW Upload 
Commands 


Monitoring and 
Analysis 


NICER SMOC 


J 
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Concept of Operations (1/3) 
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Concept of Operations (3/3) 
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• NICER in Phase C 


• Successful review progress 

— SEXTANT Baseline Review Aug 2013 

- NICER PDR Dec 2013 

- SEXTANT CDR-EPR Aug 2014 

— GCD Continuation Review Aug 2014 

- NICER IFSW Aug 2014 

- NICER CDR Sep 2014 

• Milestones 

- XFSW Build 0 (internal) Dec 2013 

- XFSW Build 1 Apr 2014 

- XFSW Build 1.1 (clean-up) May 2014 

— IFSW Build 1 integration & testing Jun 2014 
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• NICER/SEXTANT is an excellent partnership between 
science and technology (SMD & STMD) 


• SEXTANT will be a: 

- First demonstration of real-time, on-board XNAV 

- Significant historical event enabled by ISS 

- Bridge between science, technology & exploration to 
connect with the public 

• Near-term milestones 

— XFSW Build 2 Mar 2015 
— IFSW Build 2 testing 

• NICER URL 

- http : //heasarc . gsfc . nasa . gov/docs/nicer/ 
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Backup 
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SEXTANT in NICER 
Flight Software Architecture 
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